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THEORETICAL RESULTS FOR FULLY FLOODED, ELLIPTICAL 
HYDRODYNAMIC CONTACTS 

The most important practical aspect of the elastohydrody- 
namic lubrication (EHL) theory developed in the last chapter is 
the determination of the minimum film thickness within the con- 
tact. The maintenance of a fluid film of adequate magnitude is 
an essential feature for the correct operation of lubricated 
machine elements. The results for a fully flooded conjunction 
are given in this chapter and are based on the earlier work of 
Hamrock and Dowson (1977a). These results show the influence of 
contact geometry on minimum film thickness as expressed by the 
ellipticity parameter and the dimensionless speed, load, and 
materials parameters. 

In the numerical analysis the ellipticity parameter was 
varied from 1, which represents two spheres in contact or a ball 
on a plane, to 8, which represents a long ellipse and a conjunc- 
tion geometry approaching that of two cylinders or a cylinder on 
a plane in line contact. The dimensionless speed and load pa- 
rameters were varied over about two and one order of magnitude, 
respectively. Conditions equivalent to using solid materials of 


*Published as Chapter 8 in Ball Bearing Lubrication by 
Bernard J. Hamrock and Duncan Dowson, John Wiley & Sons, Inc., 
Sept. 1981. 



bronze, steel, and silicon nitride and lubricants of paraffinic 
and naphthenic mineral oils were also considered in deriving - 
numerical solutions and hence in obtaining the exponent on the 
dimensionless materials parameter. Thirty-four cases were 
considered in obtaining the fully flooded minimum- and 
central-film-thickness formulas. Experimental verification of 
these relationships is presented in Chapter 10. 

A fully flooded condition is said to exist when the inlet 
distance of the conjunction ceases to influence in any signifi- 
cant way the computed minimum film -thickness. The appropriate 
inlet dimension for the conjunction is defined as the distance 
from the center of the Hertzian contact zone to the inlet edge 
of the computing area. 

In this chapter contour plots of both pressure and film 
thickness are presented that show the essential features of 
elastohydrodynamic conjunctions. In particular the crescent- 
shaped region of minimum film thickness, with its side lobes in 
which the separation between the solids is a minimum, clearly 
emerges in the numerical solutions. These theoretical solutions 
for film thickness have all the essential features of previously 
reported experimental observations based on optical interferom- 
etry that are presented in Chapter 10. 
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8.1 Dimensionless Grouping 


The variables resulting from the isothermal EHL elliptical- 
contact theory developed in Chapter 7 are 

9 

A,B = constants used to define density of fluid, m /N 

2 

E' = effective elastic modulus, N/m 
F = normal applied load,- N 
h = film thickness, m 

Piv = asymptotic isoviscous pressure obtained from Roelands 
(1966), N/m^ 

= effective radius in x (motion) direction, m 

R^ = effective radius in y (transverse) direction, m 

u = surface velocity in x direction, m/sec 

Z-^ = viscosity-pressure index, a dimensionless constant 

2 

tiq = atmospheric viscosity, N s/m 

Dowson and Higginson (1966) found that changes in the lu- 
bricant density with pressure had little effect on the minimum 
film thickness for line or rectangular contacts. We may there- 
fore assume that the same is true for elliptical contacts. Even 
though the compressibility effect is still considered in the EHL 
theory developed in Chapter 7, the constants used to define the 
fluid in the density equation are not used in the minimum-film- 
thickness formula. Therefore the 11 variables mentioned previ- 
ously are reduced to nine, A and B being eliminated. From the 
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nine variables the following five dimensionless (groupings were 
established: 

(1) Dimensionless film thickness 

( 8 . 1 ) 


where 



( 8 . 2 ) 


The radii of curvature in equation (8.2) are shown in Figure 2.18. 
(2) Dimensionless load parameter 

W = — ^ (8.3) 


where 



(8.4) 


(3) Dimensionless speed parameter 


where 


U = 


^ 0 ^ 

E’Rx 


u 


2 


(4) Dimensionless materials parameter 


G = 


E' 

p. 

IV, as 


(8.5) 


( 8 . 6 ) 


(8.7) 


where is the asymptotic isoviscous pressure obtained 

from Roelands (1966). Blok (1965) has pointed out that the 
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asymptotic isoviscous pressure can be approximated by the 
inverse of the pressure-viscosity coefficient reported by 


Roelands (1966); thus p^^ « 1/a. 

(5) Ellipticity parameter! 


( 8 . 8 ) 


where 

a = semiaxis in the y (transverse) direction 
b = semi ax is in the x (motion) direction 

The ellipticity parameter is determined from the definition of 
the radii of curvature r^^, r^^, and >" 5 y); derivation 

is given in Chapter 3. 

The dimensionless film thickness can thus be written as a 
function of the other four parameters 

H = f(k,u,W,G) (8.9) 

In the fully flooded results presented in this chapter the in- 
fluence of the dimensionless parameters k, U, W, and G on 
minimum film thickness is studied. 

The set of dimensionless groups [h, U, W, 6, and k) is a 
useful collection of parameters for evaluating the results pre- 
sented in this and the subsequent chapter. It is also compar- 
able to the set of dimensionless parameters used in the initial 
elastohydrodynamic analysis of line contacts, and it has the 
merit that the physical significance of each term is readily 
apparent. However, a number of authors, for example, Moes 
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(1965-66) and Theyse (1966), have noted that this set of di- 
mensionless groups can be reduced by one parameter without any 
loss of generality. This approach is followed in Chapter 12, 
where the film thicknesses to be expected in each of the four 
regimes of fluid-film lubrication are represented graphically. 

8.2 Effect of Ellipticity Parameter 

The ellipticity parameter k is a function of the radii of 
curvature of the solids only rj^^, r^^, and >^ 5 y)- The 

radii of curvature in the x-direction for both solids a and b 
are used in defining the dimensionless speed and load parame- 
ters. Therefore only the radius of curvature of solid b in 
the y direction was changed in varying the ellipticity param- 
eter. In doing this the dimensionless parameters for speed U, 
load W, and materials G were held constant at the following 
values: 

U = 0.1683x10”^^, W = 0.1106x10“^, G = 4522 (8.10) 

Care was taken to ensure that the highest value of the ellipti- 
city parameter (k = 8) corresponded to conditions in the fully 
developed elastohydrodynamic region. The approach used was to 
relate the operating conditions for this near-rectangular- 
contact situation to the elastohydrodynamic region defined by 
Dowson and Whitaker (1965-66). For values of k less than 8 
the results move further into the elastic region. 
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Table 8.1 shows the computed values of minimum film thick- 
ness for 10 values of the ellipticity parameter as ob- 

tained from the EHL elliptical-contact theory developed in 
Chapter 7. From these 10 pairs of data it was possible to de- 
termine a good empirical relationship between the minimum film 
thickness and the ellipticity parameter k for the con- 

ditions considered in the computations. 

The form of relationship chosen after investigating a num- 
ber of different expressions can be written as 

1 _ ( 8 . 11 ) 

Hmin,r 

A least-squares exponential curve was fitted to the 10 pairs of 
data points to obtain values for A and B in equation (8.11). 



In addition to the least-squares fit a coefficient of determina- 
tion e was obtained. The value of e reflected the fit of 
the data to the resulting equation: unity representing a per- 

fect fit, and zero the worst possible fit. The minimum film 
thickness for a rectangular contact ^ used in equation 

(8.11) was determined by finding the value of ^ that 

gave a coefficient of determination closest to unity. The 
rectangular-contact minimum film thickness was thus deduced from 
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the present set of results for the limiting case in which the 
ellipticity parameter k approached infinity. This value of 

turned out to be 7.082x10“^ with a corresponding coef- 
ficient of determination of 0.9990, an excellent fit. The values 
of A and B in equation (8.11) obtained from the least-squares 
fit were A = 0.9966 » 1.00 and B = -0.6752 « -0.68. 

Substituting these values of A and B into (8.11) estab- 
lished the following relationship between the ellipticity param- 
eter and the minimum film thickness: 


H . 
min 


(1 


-0.68k. 

e 


( 8 . 12 ) 


where = *^min^^x '^eP'^esents the value of dimensionless 

minimum film thickness predicted by the relationship that gave 
the best least-squares fit to the numerical solutions. It is 
most significant that the computed value of A was 0.9966, or 
approximately unity, since this ensured that the minimum film 
thickness approached zero as k approached zero. Therefore, 
even though the smallest value of k used in obtaining equation 
(8.12) was unity, it appears that equation (8.12) might well be 
valid at smaller values of k since in the limiting case where 
equation k » 0 equation (8.12) satisfies the physical situation 
corresponding to a long elliptical line contact in the direction 
of motion. For the other extreme case of a large value of k a 
nominal line or rectangular-contact situation was approached, 
and the agreement with existing results was again good. From 
Dowson and Higginson (1966) the rectangular-contact minimum film 
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thickness for the dimensionless parameters given in equation 
(8.10) was 7.720x10”^, compared with 7.082x10"^ from the 
present results. The difference of only 9 percent could well be 
the result of Dowson and Higginson (1966) using an exponential 
pressure-viscosity relationship instead of the Roelands (1966) 
formula used in the present work. As pointed out in the closure 
of Hamrock and Dowson (1976b), in answer to a query from the 
late P. M. Ku, the Roelands formula for viscosity as a function 
of pressure suppresses the pressure spike somewhat and results 
in a smaller film thickness than that predicted for an exponen- 
tial viscosity-pressure relationship. 

The values of for the 10 values of the ellipticity 

parameter considered are compared with in Table 8.1. 

The percentage difference between the minimum film thickness 
obtained from the EHL elliptical-contact theory and the 

minimum film thickness from the least-squares-f it equation S 
is expressed as 



(8.13) 


Note that in Table 8.1 never exceeds +3 percent. 

Representative contour plots of dimensionless pressure are 
shown for two quite different values of the ellipticity param- 
eter k (8 and 1.25) in Figure 8.1. In these and all subsequent 
contour plots to be shown in this text, the + symbol indicates 
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the center of the Hertzian contact zone. Note that the dimen- 
sionless representation of the X and Y coordinates adopted 
in the present work caused the actual Hertzian contact ellipse 
to become a circle regardless of the value of the ellipticity 
parameter. The Hertzian contact circle is shown in each figure 
by asterisks. On each figure is a key showing the contour la- 
bels and each corresponding value of dimensionless pressure. 

The inlet region is to the left and the exit region is to the 
right in each figure. 

For k = 8 the maximum pressure is near the center of the 
contact; and even though the conditions are in the elastic re- 
gion, there is no evidence of the pressure spike that is a fea- 
ture of many theoretical solutions to the elastohydrodynamic 
problem. The pressure gradient at the exit end of the conjunc- 
tion is much larger than that in the inlet region. For 
k = 1.2 a pressure spike is visible at the exit end of the 
contact. 

Two contour plots of film tnickness for the same values of 
ellipticity parameter, 8 and 1.25, are show in Figure 8.2. 

For k = 8, Figure 8.2(a), the minimum film thickness is located 
directly behind the center of the Hertzian contact zone on the 
minor axis of the contact ellipse. For k = 1.25, Figure 
8.2(b), two minimum-film-thickness regions occur in well-defined 
side lobes that follow, and are close to, the edge of the 
Hertzian contact ellipse. These results, showing the two side 
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lobes in which tninimum-f i Im-thickness areas occur, reproduce all 
the essential features of previously reported experimental 
observations based on optical interferometry (Cameron and Gohar, 
1966). 

The variations of pressure and film thickness in the direc- 
tion of motion close to the midplane of the contact are shown in 
Figure 8.3 for three values of the ellipticity parameter. It is 
evident from Figure 8.3(a) that no pressure spike is predicted 
for k = 6, but this well-known feature of the theoretical solu- 
tion to the EHL problem is evident for k = 2.5 and k =1.25. 

An interesting feature of the film thickness profiles shown 
in Figure 8.3(b) is the reentrant form of the film shape for 
k = 1.25. This can probably be attributed to lubricant compres- 
sibility, and similar profiles were reported for line contacts 
by Dowson and Higginson (1966). The film thickness in the cen- 
tral regions of elastohydrodynamic contacts lubricated by incom- 
pressible fluids is essentially constant; but when compressibil- 
ity is considered, the local film thickness is reduced by an 
amount corresponding to the change in fluid volume with pressure. 

8.3 Influence of Speed 

If the surface velocity in the x direction u is 
changed, the dimensionless speed parameter U is modified as 
shown in equation (8.5), but the other dimensionless parameters 
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(k, W, 6 , and H) remain constant. The values at which these 
dimensionless parameters were held constant in the calculations 
performed to determine the influence of speed on film thickness 
were 

k = 6, W = 0.7371x10"^, G = 4522 (8.14) 

Values of the dimensionless speed parameter U and the 
corresponding minimum film thickness as obtained from 

the EHL elliptical-contact theory developed in Chapter 7 are 
presented in Table 8.2,. Calculations were performed for 15 
values of the dimensionless speed parameter covering nearly two 
orders of magnitude. The solutions enabled the relationship 
between minimum film thickness and the speed parameter to be 
written in the form 

^in = C 3 U ^ (8.15) 

By applying a least-squares power fit to the 15 pairs of data 
^^^i’ ^min,i^’ = !>•••. 15], the values of C 3 and C^ were 
found to be C 2 = 560 and C^ = 0.6754 « 0.68. The coefficient 
of determination e for these results was excellent at 0.9998. 
Substituting the values of C^ and C^ into equation (8.15) 
gave the values of shown in Table 8.2. The percentage 

difference between the minimum film thickness obtained from 
the EHL elliptical-contact theory and the minimum film 

thickness obtained from the least— squares fit was determined 
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from equation (8.15), and the results are given in Table 8.2. Note 
that the variation of is less than +2 percent. 

From the value obtained for and equation (8.15) the 
effect of dimensionless speed on dimensionless minimum film 
thickness can be written as 

H^n “ (8.16) 

The distributions of dimensionless pressure for two widely 

different values of the dimensionless speed parameter, 

0.8416x10“^^ and 0.5050x10“^®, are shown in Figure 8.4. 

- 1 ? 

For the lower speed of U = 0.8416x10 a pressure spike is 
evident at the exit end of the contact as shown in Figure 
8.4(a). But for the higher speed of U = 0.5050x10”^^ Figure 
8.4(b) shows that no pressure spike emerges from the calcula- 
tions. Another interesting feature of Figure 8.4 is that the 
pressure in the inlet region is higher at high speeds than at 
low speeds. 

The dimensionless film thickness contours for values of U 

of 0.8416x10“^^ and 0.5050x10“^^ are shown in Figure 8.5. 

_12 

For the lower speed of U = 0.8416x10 the minimum film 
thickness is located along a crescent-shaped region close to the 
edge of the Hertzian contact ellipse as shown in Figure 8.5(a). 

For the higher speed of U = 0.5050x10“^^ the minimun>-f ilm- 
thickness area appears in Figure 8.5(b) to be within the contact 
zone between the center and edge of the Hertzian contact ellipse. 
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The variation of pressure and film thickness in the direc- 
tion of rolling quite close to the X axis near the midplane of 
the conjunction is shown in Figure 8.6 for three values of U. 
Recall that the values of the dimensionless load, materials, and 
ellipticity parameters were held constant at the levels noted in 
equation (8.14) for the calculations in this section. In Figure 
8.6(a) the dashed line corresponds to the Hertzian pressure 
distribution. This figure shows that the pressure at any 
location in the inlet region rises as the speed increases, a 
result that is also consistent with the elastohydrodynamic 
theory for line or rectangular contacts. Furthermore, as the 
speed decreases, the height of the pressure spike decreases and 
the hydrodynamic pressures gradually approach the semielliptical 
form of the Hertzian contact stresses. Note that the location 
of the pressure spike moves downstream toward the edge of the 
Hertzian contact ellipse as the speed decreases. Similar re- 
sults to those shown in Figure 8.6(a) have been presented for 
nominal line or rectangular contacts by Dowson and Higginson 
(1966). 

A typical elastohydrodynamic film shape with an essentially 
parallel section in the central region is shown in Figure 
8.6(b). There is little sign of a reentrant region in this 
case, except perhaps at the lowest speed. There is also a con- 
siderable change in the film thickness as the dimensionless 
speed is changed, as indicated by equation (8.16). This 
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illustrates most clearly the dominant effect of the dimension- 
less speed parameter U on the minimum film thickness in 
elastohydrodynamic contacts. 

8.4 Influence of'Load 

Changes in the dimensionless load parameter W can be 
achieved while keeping the other parameters constant by changing 
only the normal applied load F in equation (8.3). The values 
at which the remaining parameters k, U, and G were held 
constant during the calculations were 

k = 6, U = 0.1683x10"^^, G = 4522 (8.17) 

Eight values of the dimensionless load parameter W cover- 
ing just over one order of magnitude and the corresponding val- 
ues of minimum film thickness as obtained from the EHL 

elliptical-contact theory developed in Chapter 7 are shown in 
Table 8.3. These eight pairs of data were used to determine an 
empirical relationship between the dimensionless load and the 
minimum film thickness, in the form 

H„. = (8.18) 

min J 

By applying a least-squares power fit to the eight pairs of data 
[(Wi, i = 1 ,..., 8 ], the values of Cg and Cg were 

found to be C^ = 2.16x10"^ and Cg = -0.0729 » -0.073. The 
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coefficient of determination e for these results was 0.9260 
and was deemed to be quite good, but it was -nevertheless the 
lowest value obtained in deriving the minimum-film-thickness 
equation (8.23). Substituting the values for Cg and Cg into 
equation (8.18) gave the values of shown in Table 8.3. 

The percentage difference between the minimum film 
thickness obtained from the EHL elliptical-contact theory 

and the minimum film thickness from the least-squares- 
fit equation was ascertained by the application of 

equation (8.13), and the results are given in Table 8.3. The 
values of never exceed +^3 percent. 

From the value of Cg and equation (8.18) the effect of 
load on minimum film thickness can be written as 


H_ 


mxn 


W 


,-0.073 


(8.19) 


This very small negative exponent is again consistent with the 
findings of earlier theoretical and experimental studies of 
elastohydrodynamic lubrication. It shows that an increase in 
load will produce a decrease in minimum film thickness but that 
the decrease is remarkably small. The elastohydrodynamic film, 
once formed, is thus quite insensitive to changes in load if all 
other factors remain constant. 

The distributions of dimensionless pressure for the two ex- 
treme values of dimensionless load parameter considered 
(W = 0.1106x10 ^ and W = 0.1290xl0~^) are shown in contour 
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form in Figure 8.7. At the lower load of W = 0.1106x10“^ 

(Figure 8.7(a)) the conditions do not promote the formation of a 
pressure spike, but this distinctive feature of many elastohydrody- 
namic solutions is evident at the higher load of W = 0.1290x10”^ 
(Figure 8.7(b)) . 

Contour plots of dimensionless film thickness for the same 
two values of the dimensionless load parameter considered in 
deriving the pressure distributions depicted in Figure 8.7 are 
shown in Figure 8.8. For the lower load of W = 0.1106x10 
the minimum film thickness shown in Figure 8.8(a) occurs di- 
rectly behind the center of the contact. For the higher load 
of W = 0.1290x10“^ the minimum film thickness is found in 
regions broadly represented by a crescent located near the edge 
of the Hertzian contact ellipse, as shown in Figure 8.8(b). 

The variation of pressure and film thickness in tne direc- 
tion of motion along a line close to the midplane of the con- 
junction is shown in Figure 8.9 for three values of dimension- 
less load parameter. The values of the dimensionless speed, 
materials, and ellipticity parameters were held fixed as de- 
scribed by equation (8.17) for all computations at various 
loads. Note from Figure 8.9(a) that the pressure at any 
location in the inlet region falls as the load increases. For 
the highest load of W = 1.106x10”^ the film thickness 
depicted in Figure 8.9(b) rises between the central region and 
the outlet restriction in the same manner as shown previously 
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in Figure 8.3(b). Again this reentrant effect is attributed to 
lubricant compressibility. Note also that at W = 0.5528x10"^ 
the film thickness is slightly smaller than at W = 1.106x10“^. 
This somewhat curious result is linked to the fact that the lo- 
cation of the minimum film thickness also changes drastically 
over this load range. At the lower load the minimum film thick- 
ness is located on the midplane of the conjunction downstream 
from the center of the contact; at the higher load it moves to 
the side lobes as described earlier. 

8.5 Influence of Material Properties 

A study of the influence of the dimensionless materials 
parameter G on minimum film thickness has to be approached 
with caution since in practice it is not possible to change the 
physical properties of the materials, and hence the value of G, 
without influencing the other dimensionless parameters consi- 
dered earlier. Equations (8.3), (8.5), and (8.7) show that as 
either the materials of the solids (as expressed in E') or 
the lubricants (as expressed in and p-^ ^^) are varied, 
not only does the materials parameter G change, but so do the 
dimensionless speed U and load W parameters. Only the el- 
lipticity parameter can be held fixed, and for all results pre- 
sented in this section k = 6 was adopted. 
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The results obtained from calculations performed for four 
values of the dimensionless materials parameter are summarized 
in Table 8.4. The general form of these results, showing how 
the minimum film thickness is a function of the dimensionless 
materials parameter, is written as 


where 



( 8 . 20 ) 


Hmin 

^ _ ^-0.68k^y0. 68^-0. 073 


( 8 . 21 ) 


In equation (8.21) the exponents have been rounded off to two 
significant figures so that any error could be absorbed in C^, 
given in equation (8.20). By applying a least-squares power fit 
to the four pairs of data, the values of and Cg were 

found to be = 3.69 and Cg = 0.487 » 0.49. The coefficient 
of determination for these results was 0.9980, which is excel- 
lent. Substituting values of and Cg into equation 

(8.20) gives the values of shown in Table 8.4. The 

percentage difference between and shown 

in Table 8.4 varies by only 2 percent at all times. Therefore 
the effect of the dimensionless materials parameter on the di- 
mensionless film thickness can be written with adequate accuracy 
as 

( 8 . 22 ) 
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8.6 Minimum-Film-Thickness Formula 


The proport ional ity expression equations (8.12), (8.16), 
(8.19), and (8.22) have established how the minimum film 
thickness varies with the ellipticity, speed, load, and 
materials parameters, respectively. This enables a composite 
minimum-f i ln>-thickness formula for a fully flooded, isothermal 
elastohydrodynamic elliptical contact to be written as 


%in = 3.63 


-0.68k, 

e ) 


(8.23) 


In equation (8.23) the constant 3.63 is different from 
Cy = 3.69 mentioned earlier to account for the rounding off 
of the materials-parameter exponent. It is sometimes more 
convenient to express the side-leakage factor in equation (8.23) 
in terms of the radius-of-curvature ratio R /r instead of 

y ^ 

the ellipticity parameter. This can be done by using equation 
(3.28), or 


k = 1. 03 ( ) 

where 

^by 


(3.28) 


(8.24) 



(8.25) 
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There is considerable merit in using equation (3.28) since 
it avoids the need to evaluate elliptic integrals of the first 
and second kinds in determining k. The minimum film thickness 
can thus be derived directly from a knowledge of the radii of 
curvature of the contacting bodies (r^^, r^^, *"by^ * 

It is interesting to compare the Hamrock and Dowson (1977a) 
elliptical-contact, minimum-film-thickness formula equation 
(8.23) with the corresponding equation generated by Dowson 
(1968) for nominal line or rectangular contacts. 

Vn.r = 2.65 ^ (8.26) 

The powers of U, 6, and W in equations (8.23) and (8.26) are 
quite similar considering the different numerical procedures on 
which they are based. It is also worth noting that the power 
of W in equation (8.23) is extremely close to the value of 
-0.074 produced by Archard and Cowking (1966) in their studies 
of elliptical contacts. 

8.7 Application of the E lastohydrodynamic Mi nimun>-Film-Thi ckness 
Formula for Line Contacts to Elliptical-Contact Situations 

The dimensionless minimun>-f i Im-thickness formula for line 
or rectangular contacts can, under certain circumstances, be 
used to estimate the film thickness in elliptical contacts with 
acceptable accuracy. Indeed, equations like equation (8.26) 
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were used for this purpose before elliptical-contact film 
thickness equations presented in this text became available. 

The problem is to determine the appropriate value of the 
load per unit length of the equivalent cylinder to be 
used in equation (8.26). The expression for employed in 
the line-contact film thickness formula is 


where F* 


W 


r 


F* 


X 

is the load per unit length of the cylinder, or 


(8.27) 


F* = I (8.28) 

The value of the applied load F on an elastohydrodynamic 
conjunction is usually known, but it is important that the basis 
for selecting the effective length of the contact should be 
specified. 

In a dry elliptical contact the pressure along the princi- 
pal axis varies according to the Hertzian formula for a semi- 
ellipse, whereas in a line contact the pressure remains constant 
in the axial direction. This at once suggests two approaches to 
the determination of an effective length of the elliptical con- 
tact; one that ensures that the maximum pressures in the ellip- 
tical and rectangular contacts are the same, and the other that 
yields the same mean pressure in the two situations. 
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8.7.1 Determination of tne Effective Lengtii of a Rectangular 
Contact Required to Generate the Same Maximum Pressure 
as That Produced in an Elliptical Contact 


The maximum Hertzian contact pressure in an elliptical con- 
tact can be written as (Chapter 3), 


= 3f 

^max 2iTab 


where 


(3.6) 


b 


/6£FR,y''^ 

yukE*/ 


(3.14) 


Likewise, the maximum Hertzian pressure in a rectangular contact 
(Harris, 1966) is given by 


= 

^max,r nX-b^ 


where 


(8.29) 


b = 
r 



(8.30) 


Thus, if the maximum pressures are assumed to be equal in the 
two cases, equations (3.6) and (8.29) can be equated to yield 


£ 



(8.31) 


When the expressions for b and b^ are introduced from 
equations (3.14) and (8.30) this becomes 


£ 



(8.32) 
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The semimajor axis a in a dry elliptical contact has been 


written in Chapter 3 as 


a = 


i 2 
6kVFR 

. irE' 


1/3 


(3.13) 


Hence 




(8.33) 


When equation (8.33) is introduced into equation (8.32), it 
is seen that the effective length of the cylinder required to 
yield the same maximum pressure as the actual elliptical contact 
is given by 


A = 1.33 


(8.34) 


For long, thin contact ellipses in which R >> R the 

y X 

magnitudes of b and b^ do not differ greatly and ^ ap- 
proaches unity, as shown in Figure 3.3. In this case A 

approaches a^ , a result that can be obtained directly from 
equatitjn (8.31). If the approximation for ^ from equation 
(3.19) is introduced, equation (8.34) can be written as 


Ji a [1.0003 + Q- 

y X 


(8.35) 
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8.7.2 Determination of Effective Length of a Rectangular Contact 
Required to Generate the Same Mean Contact Pressure as 
That Produced in an Elliptical Contact 


The mean Hertzian contact pressure in an elliptical contact 
is given by 


P 


m 


F 

irab 


(8.36) 


where b is given by equation (3.14) as before. 

Similarly the mean Hertzian contact pressure in a rectangu- 
lar contact of length i is given by 


P 


m,r 


F 

2lh 


(8.37) 


r 

where b^ is represented by equation (8.30). 

If the mean pressures are assumed to be equal in the two 
cases, equations (8.36) and (8.37) can be equated to give 


When equations (3.14) and (8.30) are introauced for b 
and b^, this becomes 


a 



(8.39) 


Once again equation (8.33) can be introduced to replace the 
first term in brackets in equation (8.39) and thus to yield 


a 


= IT 


(8.40) 
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where 


S = 1.0003 + (3.29) 

R /R 

y X 

Equation (8.40) should be used for the effective length of 
a cylinder required to yield the same mean pressure as that en- 
countered in the real elliptical contact. For large values of 
a long, thin contact ellipse is found and S ap- 
proaches unity. In this case Z approaches 1.85 a. 

For these long, thin elliptical contacts the length of the 
cylindrical contacts required to give comparable maximum and 
mean pressures to those in the real contact approach 1.33 a and 
1.85 a, respectively. 

8.8 Central-Fi Im-Thickness Formula 

In practical situations there is considerable interest in 
the central as well as the minimum film thickness in elastohy- 
drodynamic contacts. This is particularly true when traction is 
considered since the surfaces in relative motion are separated 
by a film of almost constant thickness that is well represented 
by the central value over much of the Hertzian contact zone. 

The procedure used to obtain an expression for the central film 
thickness was the same as that used in obtaining the minimum 
film thickness and will not be repeated here. The central-film- 
thickness formula obtained from the numerical results is 
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H- = 2.69 0.067^ _ g.eie 


(8.41) 


Note that the exponent on the dimensionless load parameter 
is still negative in this equation but that the numerical value 
is slightly smaller than that determined for the minimum-film- 
thickness expression. This contrasts with the numerical study 
of Ranger, et al. (1975), who found a small but positive 
exponent on W in tneir formula for the central film thickness. 

A similar expression for nominal line or rectangular con- 
junctions based on a large number of solutions to the elastohy- 
drodynamic problem for cylinders has been presented by Dowson 
and Toyoda (1979) and can be written as 


= 3.06 (8.42) 

The agreement between the predictions of equation (8.41) with 
k = «> and those of equation (8.42) is remarkably good. 


8.9 Elastohydrodynamic Lubrication of the Ball-Race Contact in a 
Ball Bearing 

High-speed, thrust-loaded, angular-contact ball bearings 
require sufficient loading to prevent gross sliding motion 
(i.e., skidding between the ball and the inner race). As dis- 
cussed in Chapter 3 (Section 3.5.2) such sliding can lead to 
surface distress and eventually to bearing destruction. The 
degree of distress depends on the adequacy of ball-race 
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lubrication, as determined by the elastohydrodynamic minimum 
lubricant film thickness relative to the roughness of the 
sliding surfaces. Harris (1971) appears to nave been the first 
to consider the effects of elastohydrodynamic lubrication in 
determining friction forces between balls and races and to 
predict the threshold thrust loading at which skidding might 
occur. The approach used by Harris (1971) is adopted here in 
association with the elastohydrodynamic lubrication results 
given elsewhere in this chapter to facilitate an analysis of the 
lubrication of ball-race contacts in ball bearings. The 
assumptions of isothermal, Newtonian behavior of the lubricant 
in the ball-race contacts used elsewhere in this chapter are 
also used here. 

The mathematical model used is based on the forces and mo- 
ments shown in Figure 8.10. The coordinate system is shown in 
Figure 8.11. Note that the ball motion at any position is 

described by the angular speeds w , , ^ and u,, about three 

Ay z 

mutually orthogonal axes through the ball center such that 

* * , , 

“jji “ Wg cos a cos e (8.43) 

* * /V 

0 )y, = 0 )g cos a sin g (8.44) 

* 

•= u)g sin a (8.45) 
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In these equations a* is the pitch angle of the ball speed 
vector relative to the x' axis, and 6* defines the yaw an- 
gle, also called the skew angle. 

The surface shear or frictional stresses occurring in the 
ball-race contacts can be defined in terms of the sliding velo- 
cities and the lubricant pressure and film thickness acting 
within these contacts. Since sliding velocities, pressures, and 
film thickness vary from point to point within the contact, the 
tangential force in any direction acting at each contact must be 
determined by integrating the shear stressesover the elliptical 
contact surface. The sliding velocities in the y' direction 
are given by 


V 

ym 


d 

2 


(1) d' 
m e 

d 


(C (i) 
mm 


u) ,)cos g “ w , sin 3 
X m z m 


(8.46) 


where m = o or i and c^ = 1 and c^ = -1. Also the 
sliding velocity in the direction of the major axis owing to the 
rotation about the y' axis is 


d 

V = V . - T w , 

''xo XI 2 y 


(8.47) 


In equation (8.46) the pitch diameter for the 
situation resulting after the forces have been applied was 
defined in Chapter 3 (equation (3.107)). Furthermore the 
contact angles 6^ and 6^ are defined by equations (3.96) 
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to (3.99), with the exception that in this section only thrust 
loads are considered, so that 6^, = 0. 

If it is assumed that the sliding velocity varies linearly 
with the distance between the surfaces, the Newtonian concept of 
viscous fluid friction reduces to 

T = n - (8.48) 

h 

The viscosity is a function of pressure as described by equation 
(5.50). Therefore the viscosity and film thickness will vary 
within the contact, and the values to be used are those obtained 
from the elastohydrodynamic analysis. The tangential forces are 
defined as 


T = 
ym 



(8.49) 


T = 
xm 



(8.50) 


Because of its orbital speed each ball must overcome a vis- 
cous drag force imposed by the lubricant within the bearing cav- 
ity. It can be assumed that the drag caused by the gaseous at- 
mosphere is insignificant; however, the lubricant viscous drag 
depends on the quantity of the lubricant dispersed in the bear- 
ing cavity. Hence the effective fluid within the cavity is a 
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gas-oil mixture that has an effective viscosity and an effective 
specific gravity. The viscous drag force acting on an orbiting 
ball can be written as 


T 

V 


pA C (d'oj Y 
V v\ e c/ 

8g 


(8.51) 


where 

Ay = drag area of ball 
Cy = drag coefficient 
g = gravitational constant 

From Figure 8.10 the following conditions of force equili- 
brium must be satisfied for steady-state operation of the bear- 
ing: 


F cos 3 
o o 


- F, cos 3. - T sin 3 + T , sin 3. - F 

i 1 xo o xi 1 C 


0 (8.52) 


F sin e - F. sin 3. + T cos 3 - T , cos 3. = 0 (8.53) 

o o 1 i xo o xi 1 

T -T^+T=0 (8.54) 

yo yi V 


In equation (8.52) is the centrifugal force acting on the 
ball and is defined by equation (3.106). Cage-to-ring and 
cage-to-ball forces have not been considered. 

The moment equilibrium conditions to be satisfied are 


J-1 J- yT? 


(a b T + a.b.T .)dX dY - . 

00 X 0 1 i XI gy 


(8.55) 
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(a b T cos B + a b. cos B. t .)dY dX= 0 
o o yo o i 1 i yi 

(8.56) 


^Vo'^yo ^ ^i^i'^yi 


- M , = 0 
gz 


where 


M , = Jti) to , 
gy c z' 


M , = Jo) 0) , 
gz c y 


(8.57) 


and J is the polar moment of inertia. As a final condition of 
equilibrium, it is recognized that for thrust-loaded bearings 


sin Bj^ + cos B^ ^ = 0 

n 


(8.58) 


Figure 8.12 indicates the relative location of the ball 
center and the race groove curvature centers for a combined ra- 
dial and thrust-loaded bearing operating at high speeds. If the 
radial displacement <s^ is set equal to zero and the re- 
sults from this change are made to Figure 8.12, trie position of 
the ball for a thrust-loaded bearing can be determined. From 
this figure the following relationships can be developed: 
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[d(f^- 0.5) 


(8.59) 


L 


2 

4 



•= 0 


cos - ^4}^ ^ ^ ^ ^t ~ ^3)^ 

- [^'(f^ - 0.5)d + = 0 (8.60) 

Equations (8.52) to (8.60) constitute a set of nine simultaneous, 
nonlinear equations that can be solved by numerical methods for 
the variables L^, L^, a.j, 6^, 6^, a*, B*, cog, and 

The test data of Shevchenko and Bolin (1957) and the analy- 
tical results of Harris (1971) are compared in Figure 8.12. 
Analytical results according to the race-control hypothesis dis- 
cussed in Chapter 3 are also indicated. Clearly the Harris 
(1971) analysis, which uses tangential forces predicted by an 
elastohydrodynamic-lubrication analysis much like that discussed 
in this section, approximates the test data more closely than 
the race-control analytical method. 

The key elements in predicting bearing performance more 
accurately are, as mentioned by Anderson (1978), better predic- 
tions of film thickness and tractive force. The elastohydrody- 
namic lubrication film thickness equations for elliptical con- 
tacts given in this chapter are thus offered to support a step 
toward the better prediction of bearing performance. Further- 
more the study of elastohydrodynamic lubrication has stimulated 
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research on lubricant rheology that is leading to the develop- 
ment of better models for calculating tractive forces. Consid- 
erable progress has thus been made in evaluating the overall 
performance of ball bearings since Harris pioneered the way in 
1971. A revised version of Harris analysis was used in develop- 
ing a computer program called SHABERTH, which incorporates ac- 
tual traction data from a disk machine. This computer program 
is used to correlate predicted and measured bearing operating 
temperatures and power loss. Since its beginning in the 
mid-1970's SHABERTH has been updated through a number of 
versions. 

An excellent discussion of the development of computer 
programs for the analysis of ball bearing performance has been 
given by Pirvics (1980). Pirvics points out that before the 
early 1960 's there were no rolling-element computer programs of 
note. For a number of years thereafter and into the mid-1970's 
efforts were made to understand and simulate particular bearing 
types as well as the nature of elastohydrodynamic lubrication. 
During this period (1965-1975) the major concern was to improve 
understanding of the basic behavior of individual contacts 
within a bearing and the way in which they affect the perfor- 
mance of a complete bearing. This led to a system analysis that 
could be integrated into the surrounding machinery and that re- 
sulted in computer programs capable of treating the time- 
transient thermomechanics of load-carrying assemblies. There- 
fore by the end of the 1970's the significance and practical 
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necessity of the system design approach for detailing tne time- 
transient, dynamic effects were recognized. 

8.10 Closure 

By using the procedures outlined in Chapter 7 the influence 
of the ellipticity parameter and the dimensionless speed, load, 
and materials parameters on minimum film thickness has been in- 
vestigated. The ellipticity parameter k was varied from 1 (a 
ball-on-plane configuration) to 8 (a configuration approaching a 
rectangular contact). The dimensionless speed parameter was 
varied over a range of nearly two orders of magnitude, and the 
dimensionless load parameter over a range of one order of magni- 
tude. Situations equivalent to using solid materials of bronze, 
steel, and silicon nitride and lubricants of paraffinic and 
naphthenic mineral oils were considered in an investigation of 
the role of the dimensionless materials parameter. Thirty-four 
cases were used to generate the following minimum- and central- 
film-thickness relationships: 

Hjjjin = 3.63 u®*68Q0.49y-0. 073(3^ _ g-0.68k) 

Hj, = 2.69 u0-67q0.53^j- 0.067^3^ _ o.61e“®*^^^) 

Contour plots have been presented that indicate in detail 
the pressure distribution and the film shape. In some solutions 


35 



the characteristic pressure spikes generated in earlier solu- 
tions to the line-contact elastohydrodynamic problem were in 
evidence. The theoretical solutions for the film thickness have 
all the essential features of previously reported experimental 
observations based on optical interferometry. 

The results presented in this chapter have permitted the 
generation of a satisfactory theoretical film thickness equation 
for elastohydrodynamic elliptic contacts operating under fully 
flooded conditions. The exponents on the various dimensionless 
parameters governing the minimum film thickness in such conjunc- 
tions are quite similar to those developed by Dowson (1968) for 
nominal line or rectangular contacts. The most dominant expo- 
nent occurs on the speed parameter, but the exponent on the load 
parameter is very small ana negative. The materials parameter 
also carried a significant exponent, although the range of this 
variable in engineering situations is limited. 

Expressions have been developed for the selection of equi- 
valent lengths of line or rectangular contacts that permit elas- 
tohydrodynamic film thickness equations for line contacts to be 
applied to elliptical contacts in some situations. The expres- 
sions developed enable comparisons to be made between line or 
rectangular and elliptical-contact situations on the basis of 
equal maximum or mean contact pressures. 

Perhaps the most significant feature of the proposed mini- 
munvfi Im-thickness formula is that it can be applied to any con- 
tacting solids that present an elliptical Hertzian contact re- 
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gion. Many machine elements - particularly rolling-element 
bearings - possess such geometry, and it is expected that the 
new minimum-film-thickness equation will find ready application 
in such fields. 

Finally a brief outline of the application of elastohydro- 
dynamic theory to the ball-race contact in a ball bearing and a 
discussion of the influence of this theory on oearing mechanics 
are presented to illustrate the usefulness of the elliptical- 
contact solutions developed in this chapter. 
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SYMBOLS 


A 

A*, B*, C*, 
D*, L*, M* , 


A 


V 


a 


a 


B 

b 

b 


C 


C 


V 



,C 


8 


c 


c 


D 

D 

D 

De 

d 

d 




constant used in equation (3.113) 
relaxation coefficients 

2 

drag area of ball, m 

semimajor axis of contact ellipse, m 

a/2m 

total conformity of bearing 
semiminor axis of contact ellipse, m 
b/2m 

dynamic load capacity, N 
drag coefficient 
constants 

19,609 N/cm^ (28,440 Ibf/in^) 

number of equal divisions of semimajor axis 

distance between race curvature centers, m 

material factor 

defined by equation (5.63) 

Deborah number 
ball diameter, m 

number of divisions in semiminor axis 
overall diameter of bearing (Figure 2.13), m 
bore diameter, m 
pitch diameter, m 

pitch diameter after dynamic effects have acted on ball, m 
inner-race diameter, m 
outer-race diameter, m 
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E 



F* 

F 

F 





f 

f 

f 

G 


b 

r 


G* 

G 

9 


modulus of elasticity, N/m 


2 


effective elastic modulus, 

2 2 

internal energy, m /s 



V 


a 


2 

a 



processing factor 


^^^min ^min^^^min^ ^ 2/2 

elliptic integral of second kind with modulus (1 - 1/k^) 

approximate elliptic integral of second kind 

dispersion exponent 

normal applied load, N 

normal applied load per unit length, N/m 

lubrication factor 

integrated normal applied load, N 

centrifugal force, N 

maximum normal applied load (at iJj = 0), N 

applied radial load, N 

applied thrust load, N 

normal applied load at angle N 

2 1/2 

elliptic integral of first kind with modulus (1 - l/k'^) 

approximate elliptic integral of first kind 

race conformity ratio 

rms surface finish of ball, m 

rms surface finish of race, m 

dimensionless materials parameter, aE 

2 

fluid shear modulus, N/m 
hardness factor 

2 

gravitational constant, m/s 
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% 

9V 

H 

H 


H 


c,s 


"n,in 


H . 
rmn,r 


u 

min,s 


"min 


H . 
mm 


h 

h. 


h. 

1 


dimensionless elasticity parameter, 
dimensionless viscosity parameter, GW^/U^ 
dimensionless film thickness, h/R^ 
dimensionless film thickness, H(W/LI)^ = 
dimensionless central film thickness, 
dimensionless central film thickness for starved 
lubrication condition 
frictional heat, N m/s 

dimensionless minimum film thickness obtained from EHL 
elliptical-contact theory 

dimensionless minimum film thickness for a rectangular 
contact 

dimensionless minimum film thickness for starved 
lubrication condition 

dimensionless central film thickness obtained from 
least-squares fit of data 

dimensionless minimum film thickness obtained from 
least-squares fit of data 

dimensionless central-film-thickness - speed parameter, 

dimensionless minimum-film-thickness - speed parameter, 

new estimate of constant in film thickness equation 
film thickness, m 
central film thickness, m 
inlet film thickness, m 
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film thickness at point of maximum pressure, where 
dp/dx = 0, m 

minimum film thickness, m 
constant, m 

diametral interference, m 

2 

ball mass moment of inertia, m N s 
integral defined by equation (3.76) 
integral defined by equation (3.75) 
function of k defined oy equation (3.8) 
mechanical equivalent of heat 

2 

polar moment of inertia, m N s 
load-deflection constant 
ellipticity parameter, a/b 
approximate ellipticity parameter 
thermal conductivity, N/s °C 
lubricant thermal conductivity, N/s °C 
fatigue life 
adjusted fatigue life 

reduced hydrodynamic lift, from equation (6.21) 
lengths defined in Figure 3.11, m 

fatigue life wnere 90 percent of bearing population will 
endure 

fatigue life where 50 percent of bearing population will 
endure 

bearing length, m 

constant used to determine width of side-leakage region 
moment, Nm 


41 



M 


M. 


M, 


m 

tn* 


m 

m 


N 


n 

n* 


n 

P 



P 

*^max 

^iv,as 


Q 

Q 


q 

q 


f 

X 


gyroscopic moment, Nm 

-Q 75 

dimensionless load-speed parameter, WU 
torque required to produce spin, N m 
mass of Dali, N s^/m 

dimensionless inlet distance at boundary between fully 
flooded and starved conditions 
dimensionless inlet distance (Figures 7.1 and 9.1) 
number of divisions of semimajor or semiminor axis 
dimensionless inlet distance boundary as obtained from 
Wedeven, et al. (1971) 
rotational speed, rpm 
number of balls 
refractive index 

constant used to determine length of outlet region 

dimensionless pressure 

dimensionless pressure difference 

diametral clearance, m 

free endplay, m 

p 

dimensionless Hertzian pressure, N/m 

o 

pressure, N/m^ 

maximum pressure within contact, 3F/2uab, N/m 

2 

isoviscous asymptotic pressure, N/m 
solution to homogeneous Reynolds equation 
thermal loading parameter 

dimensionless mass flow rate per unit width, qnQ/pQE'R^ 
reduced pressure parameter 

volume flow rate per unit width in x direction, m^/s 
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R 

R. 


■■ax' '"bx’l 
'"ay’ '"by 
'"c- *c’ ^ 
■■s’ “s’ h 


s 

s* 


s 

s 


0 



2 

volume flow rate per unit width in y direction, m /s 
curvature sum, m 

arithmetical mean deviation defined in equation (4.1), m 

operational hardness of bearing material 

effective radius in x direction, m 

effective radius in y direction, m 

race curvature radius, m 

radii of curvature, m 

cylindrical polar coordinates 

spherical polar coordinates 

defined in Figure 5.4 

geometric separation, m 

geometric separation for line contact, m 

empirical constant 

shoulder height, m 

^O^'^max 

tangential (traction) force, N 
temperature, *C 
ball surface temperature, *C 
average lubricant temperature, °C 
ball surface temperature rise, C 

(^o'PmaxV.1 

viscous drag force, N 
time, s 

auxiliary parameter 

velocity of ball-race contact, m/s 
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u 

u 

IJ 

AU 

V 

w 


w 

X 

Y 


X 


t» 


Y 


t 



y, y, y, y^ v 

z, T, 2 , 7^ 


a 


6 


B' 

r . 

Y 

Y 
Y^ 


velocity of ball center, m/s 

dimensionless speed parameter, hqU/E'R^ 

surface velocity in direction of motion, + Ujj)/2, m/s 

number of stress cycles per revolution 

sliding velocity, u^ - u^^, m/s 

surface velocity in transverse direction, m/s 

dimensionless load parameter, F/E'R^ 

surface velocity in direction of film, m/s 

dimensionless coordinate, x/R 

X 

dimensionless coordinate, y/R^ 

dimensionless grouping from equation (6.14) 

external forces, N 

constant defined by equation (3.48) 

viscosity pressure index, a dimensionless constant 

coordinate system 

pressure-viscosity coefficient of lubrication, m^/N 

radius ratio, R^/R^ 
y X 

contact angle, rad 

free or initial contact angle, rad 

iterated value of contact angle, rad 

curvature difference 

viscous dissipation, N/m^ s 

total strain rate, s~^ 

elastic strain rate, s“^ 

viscous strain rate, s~^ 
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6 

6* 




6 


r 


6 


t 


6 


X 


6 


6 


c 


e-i 

1 

^2 


n 

T) 

"^0 

■00 

e 

A 


X 


flow angle, deg 

total elastic deformation, m 

Or“l 

lubricant viscosity temperature coefficient, c 
elastic deformation due to pressure difference, m 
radial displacement, m 
axial displacement, m 
displacement at some location x, m 
approximate elastic deformation, m 
elastic deformation of rectangular area, m 
coefficient of determination 
strain in axial direction 
strain in transverse direction 
angle between ball rotational axis and bearing 
centerline (Figure 3.10) 
probability of survival 

2 

absolute viscosity at gauge pressure, N s/m 
dimensionless viscosity, n/nQ 

2 

viscosity at atmospheric pressure, N s/m 
6.31x10“^ N s/m^(0.0631 cP) 
angle used to define shoulder height 
film parameter (ratio of film thickness to composite 
surface roughness) 

equals 1 for outer-race control ana 0 for inner-race 
control 

second coefficient of viscosity 

Archard-Cowking side-leakage factor, (1+2/3 a^)~^ 

relaxation factor 
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y 



$ 


$1 

^0 

(I) 


u. 


coefficient of sliding friction 
p/ n 

Poisson's ratio 

divergence of velocity vector, (au/sx) + (av/ay) + (aw/az), s' 

? 4 

lubricant density, N s /m 
dimensionless density, p/pq 

2 4 

density at atmospheric pressure, N s^/m 

o 

normal stress, N/m 

2 

stress in axial direction, N/m 
2 

shear stress, N/m 

2 

maximum subsurface shear stress, N/m 
2 

shear stress, N/m 

2 

equivalent stress, N/m 

2 

limiting shear stress, N/m 

ratio of depth of maximum shear stress to semiminor axis of 

contact ellipse 
p„3/2 

auxiliary angle 
tnermal reduction factor 
angular location 
limiting value of V' 

absolute angular velocity of inner race, rad/s 
absolute angular velocity of outer race, rad/s 
angular velocity, rad/s 

angular velocity of ball-race contact, rad/s 
angular velocity of ball about its own center, rad/s 
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Subscripts: 

a 

b 


c 

be 

IE 

IR 

i 

K 

min 

n 

0 

PVE 

PVR 

r 

s 

Superscript: 

(~) 


angular velocity of ball around shaft center, rad/s 
ball spin rotational velocity, rad/s 

solid a 
solid b 
central 
ball center 

isoviscous-elastic regime 

isovi scous-rigid regime 

inner race 

Kapitza 

minimum 

iteration 

outer race 

piezovi scous-elastic regime 
piezovi scous-rigid regime 
for rectangular area 
for starved conditions 
coordinate system 

approximate 
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Dimensionie: 
film thickne< 
H-h/F 

A 7.08x10 

B 7.20 

C 7.40 

D 7.70 

E 8.20 

F 8.90 

G 9.80 

H 11.00 


A 4.3x10"^ 

B 4.6 

C 5.0 

D 5.5 

E 6.0 

F 6.6 

G 7.4 

H 8.2 









DimensionI 
film thickm 
H - h/R, 

A 4.25x10“ 

B 460 

C 5.00 

0 5.50 

E 6.00 

F 7.00 

G S.00 

H 10.00 


A 61.4x10" 

B 6ao 

C 63.0 

D 65.0 

E 68.0 

F 7ao 

G 78.0 

H 86.0 



Dimensionless 



Figure 8.6. - Variation qf_dimensionless pressure 
and film thickness on X axis for three values oj_ 
dimensionless speed parameter. The value of Y 
is held fixed near axial center of contact (a) 
Dimensionless pressure, (b) Dimensionless 
film thickness. 










DimensionI 
film thickni 
H • h/R, 


A 7.0x10“^ 
B 7.2 

C 7.4 

D 7.7 

E 8.1 

F 8.7 

G 9.5 

H 10.7 




Dimensionless film thickness, H • h/R„ Dimensionless pressure, 


1.8x10’^ 



Figure 8. 9. - Variation of_di mens ion less pressure 
and film thickness on X axis for three values_of 
dimensionless load parameter. The value of Y 
Is held fixed near axial center of contact (a) 
Dimensionless pressure, (b) Dimensionless 
film thickness. 






Figure 8L Id - Forces and moments acting on a bail. 



Figure 8. 11. - Coordinate system. 
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